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We show that a generalized fluctuation-dissipation relation applies upon instantaneously increas- 
ing the temperature of a deeply supercooled liquid. This has the same two-step shape of the relation 
found upon cooling the liquid, but with opposite violation, indicating an effective temperature that 
is lower than bath temperature. We show that the effective temperature exhibits some sensible 
time-dependence and that it retains its connection with the partitioned phase space visited in ag- 
ing. We underline the potential relevance of our numerical results for experimental studies of the 
fluctuation-dissipation relation in glassy systems. 



Introduction - Supercooled liquids fall out of equilib- 
rium almost inevitably. Indeed almost any liquid can be 
taken to some low-enough temperature (avoiding crys- 
tallization) where the structural relaxation processes be- 
come extremely slow and the cooling happens to be too 
quick for the system to equilibrate [HE]- In this case 
we observe that the system slowly tries to adapt to the 
new temperature and exhibits time-dependent thermo- 
dynamic properties such as energy or pressure, entering 
a regime called physical aging. In such an off-equilibrium 
situation the correlation and the response functions 
characterizing the dynamics of the liquid loose time- 
translational symmetry and the fluctuation- dissipation 
theorem (FDT) relating the two functions is not supposed 
to hold anymore [3j [4] . 

However the interesting fact about glassy systems is 
that a general fluctuation-dissipation relation (FDR), 
seems to apply when they fall out-of-equilibrium [5]- 
[9]. This FDR allows to rigorously define a time- 
dependent effective temperature that is found to reflect 
the separation of time scales occurring in the relax- 
ation. In fact the analytic solution of schematic mod- 
els [5j [lOl [11] , numerical simulations of realistic mod- 
els of glassy systems [T2HT6] and few experimental evi- 
dences [17-20 , suggest that, when the system is rapidly 
cooled (quenched) to a low-enough bath temperature 
Tb a th, the FDR may be written as follows: 

TbathdtaW) = -X{t,t')d tf C(t,t'). (1) 

In Eq.Q C(t, t') = {A(t)B(t')) is the correlation function 
of the observables A and B and x(t, t 7 ) = (A(t)}/e e ->o 
is the response function to an external perturbation e 
applied at the time t' and coupled to the system with 
an energy contribution —eB. X(t,t') on the right- 
hand side of Eq.Q is called violation factor and it al- 
lows to introduce a time-dependent effective tempera- 
ture Teff(M') = T hath /X(t,t'). Note that X = 1 at 
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short time-scales (i.e. (t — t')/t' <C 1) where the FDR 
of Eq.Q reduces to the FDT indicating that the fast 
"vibrational" dynamics of the molecules is immediately 
equilibrated at the bath temperature (T e ff = Tb at h)- Dif- 
ferently at long time-scales (i.e. (t—t r )/t' ^> 1) it is found 
that X < 1 suggesting that the slow structural dynami- 
cal rearrangements behave as if they were in equilibrium 
at a temperature T e ff = T^th/X higher than Tb at h- This 
is a well-established scenario confirmed theoretically and 
numerically in the case where the system is taken out 
of equilibrium by rapidly cooling it down to the glassy 
phase. Nevertheless slow dynamics and aging effects can 
be observed in a glassy system without necessarily cool- 
ing it. It is possible, for example, increase rapidly the 
temperature of an equilibrated supercooled liquid or a 
glass and obtain typical aging effects [2TJ [22]. In this 
case it is interesting to ask what happens to the FDR 
and the associated T e ^ since so far no study of this kind 
has been reported in the literature [34] . 

In this work we show, for the first time, that the FDR 
and the effective temperature concept can be extended 
to the case where an initially equilibrated deeply super- 
cooled liquid is taken off-equilibrium by suddenly heat- 
ing it and letting it age towards a new equilibrium state. 
We find that the FDR displays the two slopes observed 
in quenches but with X > 1 at large time-scales. This 
suggests that the slow dynamics behaves as if it were 
equilibrated at some T e ff that is lower than T& a ^. More- 
over we find that T e ff exhibits a strong dependence on 
the time elapsed from the instantaneous heating. Finally 
we show that the connection between T e ^ and the po- 
tential energy landscape found in some glassy systems is 
conserved when the same systems are quickly heated. As 
argued in the following our results suggest a new route 
to test experimentally the breakdown of the FDT and 
measure T e ^. 

Simulation - We simulate a binary mixture (80:20) 
of 1000 particles interacting via the Kob-Andersen 
Lennard- Jones (KALJ) potential [23] [24]. Equilibrium 
low-temperature configurations are produced by numer- 
ically integrating thermostatted (NVT) Newtonian dy- 
namics for times of the order of 10 8 molecular dynamics 
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(MD) steps for the lowest T investigated [35 . MD sim- 
ulations are run on high-end graphics processing units 
(GPUs) [25] that allow for a considerable speed-up of 
the simulation that typically runs ~ 20 times faster on 
GPUs than on standard CPUs. Using such a massively- 
parallel computing approach we produce from 5 x 10 3 
to more than 10 4 independent initial configurations for 
each low-temperature state from which we perform our 
numerical heating experiments. The lowest temperature 
equilibrated is T = 0.41 at density p = 1.2, where the 
mode-coupling temperature is Tmct = 0.435 [26 j. For 
each low-T configuration we follow the off-equilibrium 
dynamics after an instantaneous temperature increase (T 
up-jump) using standard Monte Carlo (MC) simulations 
in the canonical ensemble. The choice of MC dynamics 
allows to employ a zero-field algorithm that can be used 
for obtaining the integrated response function without 
actually applying any external perturbation. With this 
technique, introduced in Ref. [16] . the response is guaran- 
teed to be free of any non-linear contribution. Moreover 
this method provides the response as a function of t' (i.e. 
the instant at which the perturbation is applied) avoiding 
to run one simulation for each t r . 

Results - We study the FDR and T e ^ by building the 
parametric FD-plot reporting T^th t') versus C(t, t'). 
In this way the FDT-violation can be immediately vi- 
sualized and estimated from the parametric curve, be- 
ing X = T bath /T ef f = -T hath d X (t,t')/dC(t,t') at large 
time-scales. The variables A and £?, defining x an d C, 
are chosen to be A\^(t) = N -1 ^ rjj exp(ik • rj(t)) and 
Bk(t) = N[Ak(t) + A_k(t)], where the sum is extended 
to all N particles of the system with instantaneous posi- 
tions (ri (£),..., r/v(£)), and r]j is a bimodal random vari- 
able with zero mean. With this choice the correlation 
function coincides with the self-intermediate scattering 
function at the wave- vector k. 

Figure [I] shows that an evident FDT-violation occurs 
when Tbath is instantaneously increased from a very low 
T. Note that the FD-plot displays the two slopes already 
found in quenches, but with opposite violation. By look- 
ing at Fig. [I] it is also clear that the violation gradually 
disappears on increasing the final T. This is because the 
relaxation becomes so fast that it is not possible to see 
clearly, neither characterize, a regime where X ^ 1 in 
the FD-plot. Differently, if the final T in the up-jump 
is too close to the initial T, a very slow relaxation is 
observed which is not suitable for our study. Inset of 
Fig. [I] shows that X does not depend on the wave- vector 
k entering in the definition of the correlation and the 
response, as observed in the case of quenches [13]. More- 
over we verify that identical FD-plots are obtained if we 
increase T or if we decrease p instantaneously starting 
from initial states that are connected by the density- 
scaling equation [27] (Fig. [I]). This confirms that the 
KALJ system has well-satisfied off-equilibrium scaling 
properties in heating/expansion protocols as found for 
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FIG. 1: Response versus correlation parametric plots (FD- 
plots) for several temperature up-jumps as the final (bath) 
temperature Tbath is varied. The time t' varies between 6 x 
10 2 and t = 8 x 10 3 MC steps. A violation of the FDT 
(dashed line) characterized by X > 1 (i.e. T e ff < T) is clearly 
observed at large time-lags t — t' as evidence by the straight- 
line fits of the points departing from the FDT. It can be seen 
how, increasing the final temperature of the jump, the FDT- 
violation is gradually lost. We show also how two identical 
FD-plots can be generated by instantaneously increasing the 
volume (at fied T) or the temperature (at fixed p) in analogy 
with the findings of Ref. [15] . (Inset) FD-plot in T up-jump 
for two different wave- vectors. The fitting line of the k — 5.2 
case fits very well also the k = 7.2 case indicating that T e ff 
does not depend on the chosen wave-vector as found also in 
quenches [15] . 



cooling/compression procedures [T2| fT5] . 

In Figure [2] we show how the FDR changes as the ini- 
tial T in the up-jump is changed. Upon increasing the 
initial T the deviation from the FDT is weaker making it 
clear why low-T initial states are necessary for our study. 
However it is interesting to note that, even for an initial 
T > Tmct, the deviation from the FDT is still quite ev- 
ident. This suggests that the FDT-violation can be ob- 
served in wide range of temperatures in the supercooled 
regime above and below Tmct- 

To further characterize the time-dependence of the 
FDT-violation we choose as initial equilibrium state T = 
0.41 and as a final state T = 0.75 (p = 1.2). In this 
jump the time evolution of the system is slow enough to 
allow for a characterization of the FDR at well-separated 
aging times. Note that the time-scales and the aging 
times studied here are of the same order of those studied 
in the quenches [15]. As shown in Fig. [3] we find that 
X strongly depends on time suggesting that, as the sys- 
tem equilibrates, the deviation from the FDT decreases 
and the T e ff smoothly approaches the bath temperature 
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FIG. 2: FD-plots for T up-jumps starting with different 
initial temperatures (6 x 10 2 < t' < 8 x 10 3 MC steps). Notice 
that, as the initial temperature is increased the deviation from 
the FDT reduces. 



^bath- This dynamic behavior of the FDR gives the 
possibility of studying in detail the connection between 
T e ff and the configuration space visited out of equilib- 
rium. In the inherent state (IS) formalism[28, 29 each 
configuration of the system is associated to the config- 
uration reached by a steepest descent minimization of 
the potential energy. The phase space is partitioned in 
basins of attraction of different potential energy min- 
ima. In the IS formalism, the free-energy is written as 
F = eis - T bath s con f(eis) + /vib(eis, ^bath), where the 
configurational entropy ,s con f measures the number of dis- 
tinct minima with energy ejs while the vibrational free 
energy / v ib accounts for the contribution to the total 
free energy arising from the exploration of the sampled 
basin volume in configuration space. Such free energy 
expression is the IS analog of the TAP free energy [30 , 
where eis labels the possible TAP states. In equilibrium, 
dF/dejs = 0, providing a convenient way to calculate 
ds con f/de IS as 
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In the IS framework the aging dynamics is modeled 
as a walk in configuration space, composed by a fast 
vibration in one minimum and rare jumps among dif- 
ferent minima. This approach leads to define an off- 
equilibrium extension of the free-energy F in which the 
configurational part is weighted by the internal T of 
the system T[ n t (the analog of the effective temperature) 
while the vibrational term is assumed to be in equilib- 
rium with the bath temperature [14 . This results in 
F = eis ~ T int Sconf(eis) + /vib(eis,^bath). By analogy 
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FIG. 3: FD-plots for several T up-jumps as the time t in- 
creases. Here t' varies always from ~ t/10 to t. The deviation 
from the FDT (dashed line) gradually decreases as t increases 
suggesting that the slow degrees of freedom of system slowly 
thermalize to the bath temperature. The straight lines have 
slope — T/Tint(t) (only the intercept has been adjusted for 
fitting), where Tint is the internal temperature obtained by 
using Eq. Q Q3] (see text). 



with standard thermodynamics T- m t can be obtained by 
minimizing the free-energy with respect to ei$: 
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Substituting Eq. [2] in Eq. [3] results in a relation among 
several quantities: the eis value explored during aging, 
Tbath, the equilibrium T at which the basin of energy eis 
is typically populated (T eq ) and T int 



rri / \ rri / \ 1 + /vib(eiSj^bath) 

Ant(eis) = ^eq(eis)— — z 7 T^~ ( \T • 

1 + /vib(eis,-7eq(eis)) 



(4) 



The graphic method to evaluate T eq from the knowledge 
of the eis explored during the aging dynamics is reported 
in Fig. [4j Following Ref. [14] we determine / v ib in equi- 
librium simulations by approximating the IS minima as 
harmonic 3N-dimensional wells. The T[ n t value computed 
via Eq. Q for different aging times can be compared with 
the slope of the FDR. Results displayed in Fig. [3] show 
that, to a very good approximation, T e ff(t) = Ti nt (t), 
supporting the thermodynamic link implied in the FDR 
approach. These results suggest that, even in a T up- 
jump, T e ff is connected to the dynamics of the underlying 
IS [14j [15] . We also note that the harmonic approxima- 
tion for / v ib, that is expected to be less accurate as T is 
increased, still leads to very reasonable results. 
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FIG. 4: (Left panel): Average inherent structure energy 
in equilibrium as a function of T. (Right panel): Average 
inherent structure energy after a T up-jump as a function of 
time. The mapping between the equilibrium and the non- 
equilibrium IS energies allows to estimate the internal tem- 
perature that fits very well the FDT- violation (see Fig. [3}. 



Conclusions - We have investigated the off-equilibrium 
dynamics of a deeply supercooled atomic liquid following 
an instantaneous increase of the bath temperature. Rely- 
ing on GPU-implemented numerical simulations we have 
focused on the FDR and the effective temperature that 
can be measured out of equilibrium. We have shown that 
FDR is characterized by the typical two-step shape asso- 
ciated with the time-scale separation of glassy dynamics. 
The violation factor is found to be greater than one giv- 
ing a T e ff lower than the Tb a th temperature that is the 
opposite of what is found upon cooling a liquid in the 
glassy phase. We have identified the most suitable equi- 
librium states among which the temperature up-jump has 
to be performed to observe and characterize in detail 
the FDR. The effective temperature measured displays 
a clear time-dependence. Further analysis indicates that 
T e ff coincides to a good approximation with the inter- 
nal temperature obtained by a extended (off-equilibrium) 
thermodynamic framework based on the potential energy 
phase-space decomposition. 

Our results suggest a new route to follow in the exper- 
imental study of the FDR. In these challenging experi- 
ments one has indeed to measure independently response 
and correlation as the temperature is quickly changed. 
The temperature change has to be very fast with respect 
to the relaxation time of the system. One possibility 
suggested by our results is to try to perform an heating 
experiment among equilibrium liquid states that can be 
quite easily prepared experimentally. These experiments 
will possibly take advantage from the fact that rapidly 



heating the liquid might be much easier than cooling if 
the heating is performed by using, for example, some 
laser-based techniques [32j[33]. It could be, quite unex- 
pectedly, that the robust concepts of fluctuation dissipa- 
tion relation and effective temperature, that were con- 
ceived for characterizing the fast cooling of a liquid, are 
instead more easily accessible in a heating experiment. 

NG and FS acknowledge support from ERC-226207- 
PATCHYCOLLOIDS. 
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